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Abstract 

A novel scheme to simulate the evolution of a restricted set of observables of a quantum system 
is proposed. The set comprises the spectrum-generating algebra of the Hamiltonian. Focusing 
on the simulation of the restricted set of observables allows to drastically reduce the complexity 
of the simulation. This reduction is the result of replacing the original unitary dynamics by a 
special open-system evolution. This open-system evolution can be interpreted as a process of weak 
measurement of the distinguished observables performed on the evolving system of interest. Under 
the condition that the observables are "classical" and the Hamiltonian is moderately nonlinear, 
the open system dynamics displays a large time-scales separation between the relaxation of the 
observables and the decoherence of a generic evolving state. The time scale separation allows 
the unitary dynamics of the observables to be efficiently simulated by the open-system dynamics 
on the intermediate time-scale. The simulation employs unraveling of the corresponding master 
equations into pure state evolutions, governed by the stochastic nonlinear Schrodinger equation 
(sNLSE). The stochastic pure state evolution can be efficiently simulated using a representation 
of the state in the time-dependent basis of the generalized coherent states, associated with the 
spectrum-generating algebra. 

PACS numbers: 03.67.Mn,03.67.-a, 03.65.Ud, 03.65 Yz 
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I. INTRODUCTION 



The number of independent observables of a quantum system with the Hilbert space 
dimension N is N 2 — 1. In many-body systems, when N increases exponentially with the 
number of degrees of freedom, that large number of observables can be neither measured nor 
calculated. Only a limited number of dynamical variables is accessible to an experimentalist, 
while all the uncontrollable parameters are averaged out. This means that generically, an 
observed quantum system is characterized by a small number of the expectation values of 
accessible observables. To theoretically characterize the dynamics of a quantum system it 
is desirable: (i) to find equations of motion for this reduced set of expectation values; (ii) 
to be able to solve the associated equations of motion efficiently A computational cost of 
a direct quantum simulation scales as 0(N 5 ), 5 > 1 [1]. A simulation is defined as efficient 
if the corresponding equations of motion can be solved with a computational cost which is 
reduced substantially from that number. 

In the present study we explore the possibility of an efficient simulation of a restricted 
set of observables, using a novel paradigm for the simulation. Assuming that the set of 
experimentally accessible observables is small, it is plausible that there exist a number of 
microscopic theories, leading to the same observed dynamics. If a microscopic theory can 
be found, which leads to equations of motion that can be solved efficiently, the dynamics of 
the restricted set of observables can be efficiently simulated. More specifically, we propose 
to simulate the unitary dynamics of a quantum system by embedding it in a particular 
open system dynamics. In this dynamics the coupling to the bath is constructed to have 
a negligible impact on the evolution of the selected set of observables on the characteristic 
time scale of their unitary evolution. The key point is that the resulting open system 
dynamics can be efficiently simulated. The reduction of the computational complexity of 
the evolution, imposed by the bath, is attributed to dynamical coarse-graining, collapsing 
the system to a preselected representation which is used as the basis for the dynamical 
description. Since the bath has no observable effect by construction it should be considered 
solely as a computational tool. For that reason a term fictitious bath is used in the paper to 
refer to it. 

The quantum systems considered in the present work have finite Hilbert space dimension. 
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The dynamics is generated by the Lie-algebraic Hamiltonians: 

K = ^a i X i + J2b ij X i X j , (1) 

i ij 

where the set {Xj} of observables is closed under the commutation relations: 

K 

= fijkXk, (2) 

k=i 

i.e., it forms the spectrum-generating [2] Lie algebra [3] of the system. This algebra is 
labeled by the letter q in what follows. Lie- algebraic Hamiltonians (1) are abundant in 
molecular [4, 5], nuclear [2, 5] and condensed matter physics [2] . The basis of the algebra 
{Xj} is chosen as a distinguished set of observables, which are to be simulated efficiently 
Lie algebras considered in the present work are compact semisimple algebras [3] and the 
basis {Xj} is assumed to be ortho normal with respect to the Killing form [3]. 

The corresponding open system dynamics, which is alleged to simulate the unitary dy- 
namics of the elements of g, is governed by the following Liouville-von Neumann equation 
of motion 

—p = £p = -i |h, pj - 7 [ X i> [ X v P\ \ > ( 3 ) 

j'=i 

which has the Lindblad form [6, 7], i.e., it describes a Markovian completely positive [7] 
nonunitary evolution of the quantum system. The physical interpretation of the evolution, 
governed by Eq.(3) is the process of weak measurements [8] of the algebra of observables q, 
performed on the quantum system, evolving under the Hamiltonian (1). 

The foundation of the method is the observation that coupling to the bath induces a 
decoherence of the evolving density operator in a particular basis known as generalized 
coherent states (GCS), associated with the algebra (Section II). It is shown that if the 
Hamiltonian is linear in Xj and a certain " classicality condition" is satisfied by the Hilbert 
space representation of the algebra, the decoherence time-scale is much shorter than the 
timescale on which the effect of the bath on the elements of q is measurable, i.e., the 
relaxation time-scale. It is conjectured that this strong time separation will also hold for 
Hamiltonians bilinear in the elements of g (Section III). 

We propose to take advantage of this property of the open system dynamics for efficient 
simulation of the unitary evolution of {Xj}, using stochastic unravelings of the evolution 



[9-11] and representing the evolving stochastic pure state in the time-depending basis of 
the GCS [12, 13] (Section IV). The effect of the decoherence translates into localization 
of evolving stochastic pure state in the GCS basis, which enables efficient representation 
and simulation of the stochastic evolution. Averaging over the unravelings recovers the 
unitary dynamics of the algebra generators. It is shown that the averaging can be performed 
efficiently provided the corresponding dynamics can be efficiently measured in the lab. The 
effect of coupling to the fictitious bath is illustrated by the dynamics of a Bose-Einstein 
condensate (BEC) in a double-well trap [14, 15] modeled by the two-mode Bose-Hubbard 
Hamiltonian (Section V). It is demonstrated that the bath induces drastic localization on the 
level of a stochastic pure-state evolution, while having no observed effect on the dynamics 
of the elements of the spectrum-generating algebra of the system. 

II. EVOLUTION OF STATES 

A central theme in this section is the intimate relation between the evolution of the sub- 
algebra of observables and the dynamics of the generalized coherent states (GCS) associated 
with this subalgebra. The GCS minimize the total uncertainty with respect to the basis 
elements of the subalgebra and in addition are maximally robust to interaction with the 
bath, modeled by Eq.(3). 

A. Generalized coherent states and the total uncertainty. 

Let us assume that the subalgebra q is represented irreducibly on the system's Hilbert 
space TL. Then an arbitrary state tp G H can be represented as a superposition of the gen- 
eralized coherent states (GCS) [12, 13] \Q,ipo) with respect to the corresponding dynamical 
group G and an arbitrary state i/jq: 



where fx(il) is the group invariant measure on the coset space G/H [3] , f2 G G/H, H C G 
is the maximal stability subgroup of the reference state i/jq: 




(4) 



h\4> ) =e uKh) \4> ), heH 



(5) 
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and the GCS |0, Vo) are defined as follows: 

U(<?) |^o) = U(^) l^o) = e^U(fi) |Vo> = e^ (/l) Vo) , 9 £ G, h E H, Q E G/H, (6) 

where U(g) is a unitary transformation generated by a group element g E G. 

The group-invariant tota/ uncertainty of a state with respect to a compact semisimple 
algebra $j is defined as [12, 16]: 

W = t{**S\ = t{*5\-t {**)',■ m 

j=i j=i j=i 

The first term in the rhs of Eq.(7) is the eigenvalue of the the Casimir operator of in the 
Hilbert space representation: 

C = £x? (8) 

j'=i 

and the second term is termed the generalized purity [17] of the state with respect to g: 

p,m=e(x,v. (9) 

Let us define A mi „ as a minimal total uncertainty of a quantum state and as the eigenvalue 
of the the Casimir operator of g in the system Hilbert space. Then 

A mm < A(V>) < c H , (10) 

The total uncertainty (7) is invariant under an arbitrary unitary transformation generated 
by g. Therefore, all the GCS with respect to the subalgebra q and a reference state ipo have 
a fixed value of the total invariance. It has been proved in Ref.[16] that the minimal total 
uncertainty A min is obtained if and only if ip is a highest (or lowest) weight state of the 
representation (the Hilbert space). The value of A min is given by [16, 18] 

A min = (A, A i)<A(^)<(A,A + ^) = c w , (11) 

where A e W is the the highest weight of the representation, /i e W is the sum of the 
positive roots of q, r is the rank of g [3] and (, ) is the Euclidean scalar product in IR r . 
The corresponding CGS were termed the generalized unentangled states with respect to 
the subalgebra g [17, 18]. The maximal value of the uncertainty is obtained in states 
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termed maximally or completely entangled [17, 18] with respect to g. The maximum value 



2 



equals cu in the states having Xj ipj — for all i. Such states exist in a generic 
irreducible representation of an arbitrary compact simple algebra of observables [18]. Generic 
superpositions of the GCS have larger uncertainty and are termed generalized entangled 
states with respect to g [17, 18]. In what follows, it is assumed that the reference state ipo 
for the GCS minimize the total invariance (7). 



B. Decoherence timescales. 



The rate of purity loss in an arbitrary pure state p = \ip) {ip\ can be calculated using 
Eq.(3) as follows [19]: 

K 



d_ 

It 



Tr {p 2 } = Tr {ipp} = 2Tr ii [h, p\ p - 7 ^ [x„ [i„ p 

I 3=1 



P 




Xj,p 




X 



1> 



X, 



1> 



(12) 



i.e., the rate is proportional to the group-invariant uncertainty (7). From Eqs.(12) and (10) 
it follows that the time scale of the purity loss in a generic state is ( y r yc n )^ 1 , where is 
the eigenvalue of the Casimir, Eq. (8). On the contrary, the rate of purity loss of a GCS is 
determined by A min , Eq.(ll), which implies that GCS are rubust against the influence of 
the bath [19]. 
Assume that 



(13) 



The strong inequality (13) can be interpreted as follows. Under the action of the bath, 
modeled by Eq.(3), a generic superposition of the GCS, Eq.(4), decoheres on the fast time 
scale (7C^)~ 1 into a proper mixture of the GCS, which then follows the slow evolution on 
a time scale fixed by A min . As a consequence, the effect of the bath is to " diagonalize" the 
evolving density operator into a time dependent statistical mixture of the GCS. 

Accordingly, (7C^) _1 determines the decoherence time scale of the density operator in 
the basis of the GCS. 
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Condition (13) does not depend on the strength of coupling to the bath and therefore is 
a property of the subalgebra of observables and its Hilbert space representation. Condition 
(13) will be termed the classicality condition on the algebra of observables (see Appendix B 
for some examples). 



III. EVOLUTION OF THE OBSERVABLES 



The purpose of the present section is to show that the classicality condition (13) implies 
a large time-scales separation between the decoherence of the state and the relaxation of the 
observables comprising the spectrum-generating algebra of the system, if the Hamiltonian 
(1) is linear in the generators of the algebra. It is conjectured that the time-scales separation 
is preserved by the Hamiltonians at most bilinear in the generators. As a consequence, the 
unitary evolution of the observables on the intermediate time scale can be simulated by 
the open system dynamics, while the effect of the decoherence can be employed for efficient 
simulation of the open dynamics. 

Consider an Hamiltonian linear in the elements of the algebra g, i.e., all = in Eq. 
(1). The corresponding Heisenberg equations for the observables in g becomes: 



dt 



H, Xj 



Xj, Xj 
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K K 

k=l j,l=l 
K K 

= -ij2( ia ik)± k -iYl F*) 2 ***' ( 14 ) 
fc=i j,i=i 

where Tj fc = ifijk is a matrix element of the adjoint representation [3] of Xj. It is assumed 
without loss of generality that g is a compact simple subalgebra of observables ( in the general 
case of a semisimple algebra, the system of Eqs.(14) decouples into systems of equations for 
the simple components of the algebra). The coefficients in the r.h.s. of (14) obey 

E(^) 2 = C 2 (15) 

where C2 is the quadratic Casimir of g in the adjoint representation. Therefore 

J2( TJ )j = ( C 2h = c ad jS a (16) 
0=1 / a 



leading to 



8 , K 

—Xi = -i^2(ia ik )X k --fc ad jXi, (17) 
k=i 

which in a matrix notation reads 

= -i(A- 7 c adj )X, (18) 

where A = is defined by A ki = ia k i and X = {Xi, X 2 , ...X fc }. We define Y = 
{Y 1 ,Y 2 ,...Y fc }by 

^% = —iAYi = -iu>i%, (19) 
where uji are real since A is Hermitian. Then Y diagonalize also Eq.(18): 

—Yi = (-iA - 7c adj ) Yi = (-iui - 7c adj ) Y,, (20) 
leading to the solution of Eq.(18): 

Yi(t) = Y i (0)e- (ia;i+7Cad J ) * (21) 

and 

X i (0=X)cyY i (t). (22) 

j 

The solution (21) is obtained for an arbitrary compact simple subalgebra of the system 
observables q = su(K) C su(N) for a quantum system in a Y-dimensional Hilbert space. 
It can be generalized to a semisimple subalgebra of observables, i.e., a direct sum of simple 
subalgebras, g = ©™ =1 su(.fQ) C su(N), corresponding to a tensor-product partition of the 
system Hilbert space H = ®™ =1 Hi- In this case, Eq.(21) corresponds to local observables of 
any given subsystem. 

From Eq.(21) we see that the expectation values of observables in q oscillate on the 
timescales Ui and decay on the time scale 7c ad j. Consider an observable Y« such that 
uJi ^> 7c ad j . When the measurement of Y, in a time interval 

(uj,)- 1 <r< (7Cadj) _1 , (23) 

is performed, the nonunitary character of the evolution cannot be discovered. Therefore, 
given the time interval r any 7 with the property r <^ (7C ad j) _1 will lead to apparently 
unitary dynamics of Yj on the time interval r . 
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Next we note that since (A,/i) ^ in Eq.(ll) (a positive root has strictly positive scalar 
product with the maximal weight vector) strong inequality (13) implies |A| ^> which 
leads to the following strong inequality 

v^w > ( 24 ) 
Therefore, a time interval r exists such that 

(icnV 1 < c^r 1 < r < (TCadjT 1 . (25) 

for some i corresponding to an observable Yj in Eq.(21). The lhs of the inequality (25) is 
the decoherence rate of a generic superposition of the GCS with respect to the algebra q and 
the lhs is the decay rate of the observable Yj. This system of strong inequalities implies two 
important properties of the open system dynamics, Eqs.(14): (i) a generic superposition of 
the GCS collapses into a mixture of the GCS on a time scale much shorter than a physically 
interesting time scale of the unitary evolution of the observable; (ii) the time scale of the 
unitary evolution of the observable is much shorter than its relaxation time scale. 

It should be to emphasized that the classicality condition (13) is not sufficient to imply 
properties (i) and (ii) in a generic case of nonlinear Hamiltonians (1). Generic nonlinearity 
in the Hamiltonian is expected to introduce faster time-scales in the relaxation of the ob- 
servables. The crucial question is whether these timescales are as short as the decoherence 
timescale (7c w ) _1 . Since the decoherence time-scale depends not only on the algebra but 
also on its Hilbert space representation, it seems that in order to introduce as fast time-scales 
into the relaxation, the nonlinearity should effectively couple most of the Hilbert-Schmidt 
basis of the system operators on the physically interesting time-scale of the unitary evolu- 
tion of the algebra elements. It is conjectured that a Hamiltonians at most bilinear in the 
elements of the spectrum-generating algebra does not possess such strong coupling property. 
A closer investigation of this important question is left to future research. 

IV. EFFICIENT SIMULATION OF THE EVOLUTION OF THE SPECTRUM- 
GENERATING ALGEBRA OF OBSERVABLES. 

Efficient simulation is defined as a simulation based on a numerical solution of the first 
order differential equations for a number of dynamical variables which is much smaller than 
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the Hilbert space dimension of the system. The number of dynamical variables m cannot 
be smaller than the number of observables to be simulated, which equals the dimension K 
of the spectrum-generating algebra q. If there is a large gap between the dimension of the 
algebra and the Hilbert space dimension K = dim{g} <C dim{7i} = N the simulation based 
on the the number of variables K < m <C iV is considered efficient. 

The principle of efficient simulation of the observables, forming the spectrum-generating 
algebra g of the Hamiltonian (1) is based upon 

• Simulating the unitary evolution of the observables by the fictitious open system dy- 
namics, governed by the Liouville-von Neumann Eq.(3); 

• Unraveling the Liouville-von Neumann Eq.(3) into pure state evolutions, governed by 
the stochastic nonlinear Schrodinger equation (sNLSE) (see below); 

• Efficient simulation of the stochastic nonlinear pure state dynamics, using expansion of 
the state in a time-dependent basis of the generalized coherent states (GCS), associated 
with the spectrum-generating algebra g. 

In the previous section we have discussed the first of the listed items. The other two items 
focus on the principles of efficient simulation of the open-system evolution. 

Solving directly the Liouville-von Neumann master Eq. (3) is more difficult than the origi- 
nal problem. A reduction in complexity is based on the equivalence between the Liouville-von 
Neumann equation and the stochastic nonlinear Schrodinger equation (sNLSE)[9-ll]: 

d\rl>) = j-^-7^(x,-(x,)J dt + J2(^i-(^i)^jd^m, (26) 
where the Wiener fluctuation terms d£i satisfy 

< d£i >= 0, didij = 2-fdt. (27) 

To demonstrate the equivalence, Eq.(26) can be cast into the evolution of the projector 

= IV>> M 

dP* = (-i [H, FV] - 7 J2 *V] ] ) dt + E { " (^*) J d ^ P ^ } • ( 28 ) 

Averaging Eq.(28) over the noise recoveres the original Liouville-von Neumann equation 
(3). Therefore, the problem of efficient simulation of the Liouville-von Neumann dynamics 
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is transformed to the problem of efficient simulation of the nonlinear stochastic dynamics, 
governed by sNLSE Eq. (26). 

The simulation of the pure state evolution according to the sNLSE(26) is based on an 
expansion of the evolving state in the time- dependent basis of the GCS, Eq.(4). In the case 
of a finite Hilbert space an arbitrary state can be represented as a superposition of M < N 
GCS: 

M 

|V) = 5>|^,A>, (29) 
i=i 

where f2j is an element of the coset space G/H, G is the dynamical group of the system 
generated by g, H is the maximal stability subgroup, corresponding to the reference state 
| A) and A is the highest weight of the Hilbert space representation of the algebra. The coset 
space G/H has natural symplectic structure [13] and can be considered as a phase space of 
the quantum system, corresponding to q. Accordingly, fi, is a point in the phase space. The 
total number of variables defining (up to an overall phase) the state ip (29) equals M times 
the dimension of the phase space G/H plus the number M of amplitudes q. The dimension 
of G/H depends on the properties of the Hilbert space representation of the algebra, but is 
always strictly less then the dimension of the algebra K [13]. Therefore, the number m of 
real parameters, characterizing the state ip (29) satisfies the following inequality 

m<M(K + 2). (30) 

It follows that the necessary condition for efficient simulation of the dynamics is that 1 < 
M <^ N in the physically relevant time interval. 

It is assumed that initial state of the system is a GCS, corresponding to M — 1 in the 
expansion (29). If one omits the nonlinear and stochastic terms in the Eq.(26) it becomes 
a regular Schroedinger equation, governing the unitary evolution of the state. If the Hamil- 
tonian in Eq.(26) is linear in the elements of q, the initial GCS evolves into a GCS by the 
definition, Eq.(6). Restoring the nonlinear and stochastic terms to Eq.(26) breaks the uni- 
tarity of the evolution but a GCS still evolves into a GCS under the full equation, Ref.[20]. 
Therefore a GCS solves the sNLSE(26), driven by a linear Hamiltonian. In Ref.[20] it is 
proved that a CGS is a globally stable solution in that case, i.e., an arbitrary initial state 
evolves asymptotically into a GCS. 
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Adding bilinear terms to the Hamiltonian (1) breaks the invariance of the subalgebra g 
under the action of the Hamiltonian and, as a consequence, an initial GCS evolves into a 
superposition of a number M > 1 of the GCS (29) in the corresponding unitary evolution. 
If the number of terms M becomes large, M = O(N), the unitary evolution can no longer 
be simulated efficiently The nonlinear and stochastic terms (representing the effect of the 
fictitious bath) in Eq.(26) is expected to decrease the effective number M of terms in the 
expansion (29) of the evolving state. This effect will be termed localization. The natural 
measure of the localization is the total uncertainty of the evolving state with respect the 
spectrum-generating algebra g or, equivalently, the generalized purity of the state with 
respect to g [21]. 

The localizing effect of the bath is proved in the Ref.[20]. Heuristically, it can be under- 
stood as follows. If each sum in the sNLSE (26) is replaced by a single contribution of a 
given operator X the uncertainty of the evolving state with respect to X is strictly decreasing 
under the action of the bath, unless the state is an eigenstate of X, in which case it vanishes 
[9-11]. Therefore, the effect of the bath is to bring an arbitrary state into an eigenstate of 
X. In our case, the observables X, are noncommuting and cannot be diagonalized simul- 
taneously. Therefore, it is expected that the effect of the bath in this case will be to take 
an arbitrary state to the state which minimizes the total uncertainty with respect to the 
elements of the algebra, i.e., to a GCS. 

The characteristic time scale of the localization is the decoherence time scale (70^) -1 . 
If the classicality condition (13) holds and the nonlinearity of the Hamiltonian is moderate 
(Cf. the end of Section III), the localization is effective on a time interval much shorter 
than the relaxation of the observables in g. As a consequence, the unitary dynamics of 
these observables can be obtained by (i) simulating the nonlinear stochastic evolution of 
the localized pure states, (ii) calculating the expectation values of the observables in each 
stochastic unraveling and (iii) averaging over the stochastic realizations. 

Calculating the expectation values and averaging (steps (ii) and (iii) above) are not part 
of the definition of efficient simulation, and therefore should be considered separately. Even 
if the step (i) can be performed efficiently according to the definition, it is left to show that 
the complexity of performing steps (ii) and (iii), measured, for example, by a number of 
elementary computer operations, scale substantially less than the size of the Hilbert space 
dimension. 
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To calculate the expectation value of an observable in a state represented by the GCS 
expansion (29) one has to calculate M(M + l)/2 matrix elements of the operator between 
the GCS. Each matrix element for an operator X; G g can be calculated group-theoretically 
[13, 22], i.e., independently on the Hilbert space representation. Therefore, if M -C N the 
computation of the expectation values of the elements of g can be performed efficiently 

Complexity of the step (iii) is measured by the number of stochastic realizations necessary 
to obtain the expectation values of the observables to a prescribed accuracy. In Appendix C 
it is shown that the averaging can be performed efficiently, provided the expectation values 
of the elements of g can be measured efficiently. More precisely 

n st {e) < n ex (e)dim{g}, (31) 

where n st (e) is the number of stochastic realizations, necessary to obtain the expectation 
value of each observable Xj e g to an absolute accuracy e, n ex (e) is the number of experi- 
mental runs, necessary to obtain the expectation value of each Xj to the absolute accuracy e 
and dimjg} is the dimension of the subalgebra of observables which is assumed to be a small 
number. It follows that the averaging can be done efficiently, provided the measurement 
can be performed efficiently, which is assumed. In addition, it is important to emphasize 
that it is not necessary to converge the averaging process in order to obtain a meaningful 
information: even a single "trajectory" bears important information. 

We shall finally focus on the step (i) of simulating the nonlinear stochastic evolution of 
the localized pure states. The localization means that the number of GCS terms M in the 
expansion (29) is much smaller than the Hilbert space dimension N and, by virtue of the 
inequality (30), the number m of parameters that characterize the evolving state is much 
smaller than N. 

The details of the derivation of equations of motion for the parameters will be given 
elsewhere [23]. Here we point out the main ingredients of the derivation. We put the sNLSE 
(26) in the equivalent exponential form 

K / N 2 



|V > +\d^j > = exp 1-iKdt - (Xi - (Xi) J dt + J2 ( X * " ( X *) J d & 




|V > 
(32) 

IV >, 
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using the fact that the infinitesimal transformations commute to the leading order. 
The transformation 



\ip' >= e 



-mdt 



W > 



(33) 



is a unitary evolution, corresponding to the Schroedinger equation. The first order differen- 
tial equation of motions of parameters of the representation (29) under this unitary evolution 
can be derived variationally [24], using (29) as a variational ansatz. Therefore, the unitary 
evolution can be simulated efficiently, provided the number of terms in the expansion (29) 
is small. 

Consider the second, nonunitary transformation 



where the starred equality follows from the fact that the Casimir operator ^2 i=1 X? act 
as identity on an arbitrary state tp', and the double-starred equality follows from the fact 
that a not-necessarily-unitary transformation generated by an element of the algebra maps 
a GCS to a GCS modulo a complex phase [13]. This transformation can be performed 
group-theoretically [13], i.e., efficiently. 

The unitary evolution, Eq.(33), generated by the nonlinear Hamiltonian (1), will lead to 
derealization of the evolving state. The nonunitary evolution, Eq.(34), — to localization. 
At sufficiently strong localization the number of terms M necessary to converge the solution 
of the sNLSE (26) on a fixed time interval will be much smaller, than in the corresponding 
unitary evolution, and efficient simulation of sNLSE (26) will become feasible. 

The next section takes up an example of a two-mode Bose-Hubbard model of a Bose- 
Einstein condensate in a double-well trap to illustrate the localizing properties of the ficti- 
tious bath. 




M 



M 




(34) 



i=i 



i=i 
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V. EXAMPLE: TWO-MODE BOSE-HUBBARD MODEL. 



A common model for an ultracold gas of bosonic atoms in a one-dimensional periodic 
optical lattice is described by the Bose-Hubbard Hamiltonian: 

H = - A + + \ X>^) 2 ' (35) 

i i 

where A is the nearest neighbors hopping rate and U is the strength of the on-site interactions 
between particles. In the simplest case of a two-sites lattice model, which has been realized 
experimentally by confining a condensate in a double- well trap [14, 15], the Hamiltonian 
(35) reduces to 

H = -A(a t 1 a 2 + aj&i) + ^ ((a^) 2 + (a 2 a 2 ) 2 ) , (36) 
where A is the tunneling rate. Eq. (36) can be transformed [25] to the su(2) set of operators 

= -(aja2 + a!,ai) 

J y = Yi^ 2 ~ ^ 

J z = — — &2&2) 

leading to the following Lie-algebraic form 

H = -luJ x + m 2 z , (38) 

where u = 2 A. The Hilbert space of the system of iV bosons in this model corresponds 
to the j = N/2 irreducible representation of the su(2) algebra. We seek to simulate the 
evolution of the operators (37), driven by the Hamiltonian (38), where the initial state of 
the system is a GCS with respect to the su(2), the spin-coherent state [13, 26, 27] . More 
specifically, the initial state is chosen as 

W0)> = |-i>, (39) 

which corresponds to the state of the condensate, localized in a single well. 

The dynamics driven by the weak measurement of the operators (37) on the evolving 
condensate is described by the Liouville-von Neumann equation of the form (3): 

d 2 

-p=-^[H,p]- 7 ^[J,,[J„p]]. (40) 

i=0 

15 
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FIG. 1: The purity and expectation values of observables as a function of time. An initial GCS, 
Eq.(39), undergoes (i) unitary, 7 = (solid lines); (ii) nonunitary, 7 = 0.05/j (dashed lines), 
evolution according to the Liouville Eq. (40). The values of parameters chosen for the numerical 
solution are uj = 15 and U = to/2j. The observed dynamics of the expectation values of J x /j, jy/j 
and J z /j is negligibly affected by the bath while the generalized purity -P SU (2)[V'] °f t ne stochastic 
unraveling of the nonunitary evolution is larger by the factor of 15 than the minimal purity of the 
unitarily evolving state. 

The classicality condition (13) for the 2j + 1 = N + 1-dimensional representation of 
the 5u(2), corresponding to N atoms in the trap, translates into the N ^> 1 condition 
(Appendix B). Therefore, for sufficiently large numbers of atoms in the trap the classicality 
is satisfied and a sufficiently weak measurement of the operators J x , J y and J z is expected to 
induce strong decoherence in the spin-coherent state basis, but leaving the dynamics of the 
operators practically unperturbed. As a consequence, the generalized purity of a stochastic 
unraveling of the Eq.(40) P su {2)[^} = J2i * s expected to remain close to unity, which 

enables efficient simulation of the corresponding dynamics. 

Fig. (1) displays the evolution of the expectation values of the operators J x /j, J y /j and 
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FIG. 2: Generalized purity averaged over a small number (2-10) of stochastic unravelings of the 
Liouville-von Neumann Eq. (40). Initial state and parameters of the equation are as in the Fig. 
(1). Purity is plotted for j = 2,4,8, 16,32,64, 128. The insert shows the generalized purity as a 
function of At larger j the value of the averaged purity is apparently consistent with the 
estimate 1 — jf(ui a U b, y c ), with f(ui a U b j c ) = 3, corresponding to M = 3 number of the GCS terms 
in the expansion of the solution. 

J z /j in the unitary evolution 7 = and in the nonunitary case 7 = 0.05/j for N = 2j = 128 
particles in the condensate. The hopping rate uj = 15 and the strength of the on-site 
interaction is U = ui/2j. It can be seen that the evolution is negligibly perturbed by the 
bath for the chosen strength of the coupling 7. We also plot the generalized purity of the 
unitarily evolving state and of a random stochastic unraveling of the nonunitary evolution. 
The generalized purity in the unitary case decreases to the value of about 0.06, which 
corresponds (Appendix A) to the number of configurations M = 0.75(2j + 1) « 100 = O(N) 
in the GCS expansion of the solution. On the other hand, the generalized purity in the 
stochastic unraveling is about 0.9 — 0.95 which corresponds to a drastic reduction of the 
number of configurations to M — 0.04(2j + 1) w 5 N. 
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An interesting feature of the stochastic evolution displayed on Fig.(l) (and observed 
in other numerical simulations, see Fig. (2)) is that apparently, the generalized purity ap- 
proaches a constant value on average. Since the generalized purity is a measure of localization 
of the state on the corresponding phase space (which is the Bloch sphere for the su(2) al- 
gebra [13, 26, 27]) such behavior is suggestive of a soliton-like solution of the sNLSE (26). 
Investigation of existence and properties of these soliton-like solutions seems to be an inter- 
esting topic for future research. For the time being let us assume that the stationary (on 
average) value P of the generalized purity as displayed on Fig.(l) is an analytical function 
of 1/j (see Fig. 2 for some evidence). Then 



to the lowest order in 1/j, where / is an unknown function of the dimensionless argument 
uo a U b Y and a + b + c = 0. Using the estimate (Appendix A) for the number of configurations 
we obtain 



i.e., the number of configurations in the expansion of the stochastic unraveling does not de- 
pend on j. Numerical evidence implies that generally f(uj a U b/ y c ) 7^ 1- For example, the value 
of f(u a U b Y) deduced from the Fig. 2 is three. This implies, that asymptotically, as j — > oo, 
the dynamics of the single-particle observables of the two-modes Bose-Hubbard model can 
be reproduced not by an averaging over stochastic GCS evolutions (stochastic mean-field 
solutions), but rather by an averaging over the stochastic evolutions of superpositions of a 
constant small number M > 1 of GCS. 

The current observations were also found in different parametric regimes of the Bose- 
Hubbard model. Similar behavior has been observed in the study of the su(2)-Hamiltonians, 
such as the Lipkin-Meshkov-Glick model [28] of a system of interacting fermions. 

VI. DISCUSSION AND OPEN QUESTIONS 

The novel strategy for efficient simulation of a unitary evolution of a restricted set of 
observables has been outlined (Cf. Fig.3). The restricted set comprises the spectrum- 
generating Lie-algebra of the system. The unitary evolution of the observables is simulated 



P = l--f(uj a U b 1 c ), 



(41) 



M = (2j + 1)(1 - VP) = f(u a U b -f c ), 



(42) 
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by a particular open-system dynamics, corresponding to the process of weak measurement 
of the observables, performed on the evolving quantum system. The scheme is based on a 
large time-scales separation between the decoherence of the evolving state in the basis of the 
generalized coherent states (GCS), associated with the algebra, and the relaxation of the 
elements of the algebra. The necessary condition for the successful implementation of the 
algorithm is the " classicality" condition on the spectrum-generating algebra and its Hilbert 
space representation. This time-scale separation is proved for linear Hamiltonians and it is 
conjectured that the time-scales separation is preserved if the Hamiltonian is at most bilinear 
in the elements of the algebra. Numerical evidence obtained in the su(2) case supports the 
conjecture. The "classicality" condition excludes efficient simulation of certain subalgebras 
of observables (Appendix B). For example, the unitary dynamics of local observables of a 
composite system of qubits cannot be simulated with higher efficiency by the open-system 
evolution. 

The fast decoherence reduces the computational complexity of the evolution, while the 
slow relaxation leaves the dynamics of the restricted set of observables practically unaf- 
fected on physically interesting time-scales. The effect of a fictitious coupling to a bath can 
be viewed as a dynamically induced coarse-graining of the evolving state in the phase-space, 
associated with the spectrum-generating Lie-algebra. The fine structure of the evolving 
state, irrelevant for the expectation values of the "smooth" observables, is rubbed out by 
the decoherence, thereby reducing the computational complexity of the evolution. This 
coarse-graining can be seen as a generalization of the process of conversion of quantum 
correlations (entanglement) to classical correlations under the action of local dephasing en- 
vironments [29] . The reduction of the computational complexity is realized by simulating the 
stochastic nonlinear Schrodinger equation (sNLSE), governing the stochastic unraveling of 
the nonunitary evolution. The GCS are globally stable solutions of the sNLSE, correspond- 
ing to an Hamiltonian, linear in the algebra elements [20]. Numerical evidence obtained 
in the su(2) case suggests that Hamiltonians bilinear in the generators asymptotically lead 
to soliton-like stable localized solutions of the corresponding sNLSE. Averaging over the 
stochastic realizations of the open-system evolution recovers the unitary dynamics of the 
restricted set of observables. It is proved that the averaging can be performed efficiently, 
provided that the observables can be efficiently measured in an experiment. 

The fictitious bath is fine-tuned — it corresponds to a process of weak measurement of the 
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Effective Hamiltonian 

H=IaiXi +Ibij XiXj 



Selected observables <Xj> 



GCS X^ |Q,\|/ > 
representation 



Fictitious Bath 

-yI[Xj,[Xj, P ]] 



Simulating the bath by 
the stochastic non-linear 
Schrodinger equation for \\f k 



Averaging 

<Xj> u =<Xj> st = i^KvkiXjtyk) 



FIG. 3: A schematic flow chart of the proposed approach to simulate dynamics of the operators 
Xj of the spectrum-generating algebra of the system. The unitary evolution of the observables is 
simulated by the open-system evolution, modeling weak measurement of the evolving observables. 
The open-system dynamics is unraveled into stochastic pure-state evolutions, efficiently simulated 
using expansion of the pure state in the GCS. Efficient averaging over n s t realizations obtains the 
expectation values, corresponding to the unitary evolution /Xj 



si 



orthonormal basis set of the operators, performed with equal rates and strengths. This fine- 
tuned bath is constructed as a computational tool. On the other hand, if the fine-tuning 
condition is dropped, the resulting open-system dynamics can represent a real physical 
situation, where the linear part of the Hamiltonian is perturbed by the time-dependent 5- 
correlated noise [30]. In that case the density operator of the system will follow an open 
evolution, corresponding to a process of weak measurement of the algebra elements, per- 
formed with generally different rates [30]. It is expected, that if the noise is sufficiently 
weak, the constant part of the Hamiltonian will induce fast (on the relaxation time-scale) 
rotation in the Hilbert-Schmidt operator space, which effectively will average out the dif- 
ference between the contributions of various measurements. Therefore, this real bath is 
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expected to induce the same type of localization as the fine-tuned fictitious bath. Numerical 
evidence obtained in the su(2) case supports this conjecture [23]. Restricting the measure- 
ments to the algebra elements, the experimentalist will not observe the effect of the bath 
if the noise is sufficiently small, while measuring the higher order correlations will reveal 
the nonunitary character of the evolution. Generally, it is expected that the resulting open- 
system dynamics can be simulated with higher efficiency than the corresponding unitary 
dynamics, provided the " classicality" condition holds. 
The main open questions to be addressed are: 

• Investigation of the effect of nonlinear terms in the Hamiltonian (1) on the relaxation 
time-scales of the observables in the spectrum-generating algebra in the corresponding 
fictitious open-system dynamics, Eq.(3). 

• Development of an efficient and convergent algorithm for simulating the evolution of 
a state in the GCS basis representation. 

• A rigorous proof of the conjecture, that for Lie-algebraic Hamiltonians, bilinear in the 
generators, Eq.(l), the generalized purity of the stochastic unraveling of the corre- 
sponding open system evolution, Eq.(3), is stationary on average and the stationary 
purity approaches a limiting value independent on the dimension of the Hilbert space 
if the classicality condition becomes stronger. 
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APPENDIX A: RELATION OF THE GENERALIZED PURITY TO THE NUM- 
BER OF CONFIGURATIONS IN THE GCS EXPANSION OF THE STATE. su(2) 
CASE. 



The phase space of a quantum system, associated with the su(2) spectrum generating 
algebra is a two-dimensional sphere [13, 26, 27], usually called a Bloch sphere. The localiza- 
tion of a state of the system in the phase space means localization of its P-distribution 
[13, 27] about a point in the phase space. Without loss of generality it can be assumed that 
the state is localized about the origin: an appropriate unitary transformation, generated by 
the su(2), maps a state localized about an arbitrary point to the state, localized about the 
origin, leaving both the generalized purity and the number of the GCS in the expansion 
invariant. For definiteness let us assume that the P-distribution has a finite support area S 
of radius a about an origin on the phase space. Using the expression for the resolution of 

identity in terms of the GCS [13, 27] |t) 

2j + 1 f d 2 r 



I = 



\r) (r\ 



7T J (l + |r| 2 ) 2 

the number of the GCS in the expansion of the state can be estimated as follows 



a | ■ 



d\r\ 



1 + W 



2\2 



(2j + 1) 
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2 


1 + 


a 


2 



(Al) 



(A2) 



To calculate the generalized purity we must calculate the expectation values of J x , J y and 
J 2 . Given the P- representation of the state, the expectation value of an observable X can 
be calculated using its Q-representation: 



X = 



2j + l 



IT 



d 2 T 



where Q±(t) — It 



X 



T 



(i + M 2 ) 2 

We have [13, 27] 

.r + r* 



P(t)Q*(t), 



(A3) 



Qh = j 



1 + |r| 2 ' 
t-t* 
i(l + \r\ 2 ) 
M 2 -l 



(A4) 



1 + |r| 2 ' 

Assuming that P(r) is symmetric about the origin (r = 0), we see that the expectation 
values of J x and J„ vanish and 



2j + l 



TV 



d 2 T 



;i + |r| 2 ) 2 
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P(r)j 



\t\ 2 -1 
1 + |t| 2 ' 



(A5) 



We assume that 

{», Irl < lad 
(A6) 
0, \t\ > \a\. 

The distribution (A6) as it stands does not correspond to a pure state. Nonetheless, it can 
be understood as a coarse grained version of a localized pure state, useful for calculation 
of the expectations of J x , J y and J z and the generalized purity P su (2)[ip], Eq.(AlO). In fact, 
Eq.(A4) gives the characteristic scale of unity for the change of the Q representation in the 
integral (A3). On the other hand, the resolution of identity (Al) implies the characteristic 
scale of the fine structure of the P-distribution (the width of the overlap of two coherent 
states) of the order of (1 + \t\ 2 )/^/j. Therefore, as long as (1 + \a\ 2 )/y/j <C 1 in Eq.(A6) 
the coarse grained distribution can be used for calculation of the generalized purity. As can 
be seen below, Eq.(AlO), for j > 1 the coarse grained description is valid for calculation of 
the generalized purity asymptotically as 1/j. 

For a particular form of the distribution (A6), Eq.(A5) simplifies to 

■l«l 2 
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(A7) 



■'-"■' { ' 2 '- l) l ^^ = J - p3{23 + 1) V-0TJ^) 

The number p in Eq.(A6) can be found from the normalization condition: 

/a 2j + l f d 2 r 7-i / \ , x f H2 d\r\ 2 

= p(2j + 1) TTW- (A8) 

from which p — (1 + |a| 2 )/(|a;| 2 (2j + 1)). Inserting this expression into Eq.(A7), we obtain 



(A9) 



1 + a 



Therefore, 



P BU(2)[ 



and 



^m'-^-y-yw) <aio) 



\a\ 2 = ^===-l (All) 
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Inserting the latter expression into Eq.(A2), we obtain for the number of GCS in the expan- 
sion: 

MM = (2j + 1) (l - y/P^]) ■ (A12) 
As argued after the Eq.(A6) expressions (A10) and (A12) are valid for P su (2)[4>] 3> l/j. 

APPENDIX B: CLASSICALITY CONDITION: (I) SUBALGEBRA su(n) OF SIN- 
GLE PARTICLES OBSERVABLES OF THE n-MODES BEC IN AN OPTICAL LAT- 
TICE; (II) SUBALGEBRA OF LOCAL OBSERVABLES OF A SYSTEM OF n d- 
LEVEL SYSTEM. 

1. BEC 

The spectrum-generating algebra of the Bose-Hubbard model for the n-modes BEC in 
optical lattice is su(n) subulgebra of the single particles observables [31, 32]. It is shown the 
the classicality condition (13) is satisfied in this case, provided the number of atoms N in 
the condensate complies with 

A>n. (Bl) 

The Hilbert space of the condensate is a totally symmetric irreducible representation of the 
su(n) [N] [5] and the value of the Casimir in this representation is [5] 

n — 1 

c n = ^—N(N + n). (B2) 
An 

The total uncertainty in the GCS by [16, 18] 

n — 1 1 

A mm = c H - (A N \A N ) = c n - — A 2 = -N(n - 1), (B3) 

where we have used the known expression [5] for the norm of the maximal weight vector [3] 
A^r in the totally symmetric irreducible representation of the su(n) [N]. The value of the 
Casimir in the adjoint representation is [5] 

Cadj = n. (B4) 
Thus Eq.(13) holds if and only if Eq.(Bl) holds. Moreover, 



;,/ X A(A + n), (B5) 



c adj V 2ri 2 
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which implies Eg. (23), provided Eq.(Bl) holds. 

Therefore, using the stochastic NLS Eq.(26), propagation can be advantageous for calcu- 
lation of the single particles observables, provided the on-site interaction preserves the time 
scale separation in Eg. (23). 



2. Local observables 



Let g be a subalgebra of local observables on the composite Hilbert space. For simplicity, 
let us consider n d level systems in the Hilbert space H = ®™ =1 Hi and a subalgebra of 
local observables q = ©™ =1 su(L) C (Bf =1 su(d) C su(d n ). Since the minimum of the total 
uncertainty (7) for a local subalgebra is obtained in a product state ip pro d = ©ILi^i; where 
each ipi is a GCS with respect to the local subalgebra su(L) it follows that 

n 

A min = A[ip prod ] = c H - P 5 [ip prod ] = ^2 ( c n t ~ P S u{L)[A\) = n (c Hd - P su (l)[GCS]) 

i=i 

= nA d , min , (B6) 

where Ti d is the Hilbert space of a rf-level subsystem and A dymin is the minimal total uncer- 
tainty of a state of any subsystem with respect to the subsystem subalgebra su(L) . Therefore, 
the condition (13) is equivalent to 

A m i n A d ^ m l n /T3 r 7\ 

= < 1, (B^J 

c n c Hd 

i.e., holds if and only if the local subalgebras su(L) of the subsystems operators comply with 
the classicality condition. For example in the composite system of two-level system the only 
subalgebra of local observables is the local subalgebra $j = ®" =1 su(2). The eigenvalue of 
the local Casimir equals (1/2) (1/2 + 1) = 3/4 and the generalized purity with respect to a 
su(2) algebra of each two- level system is 1/4. Therefore the minimal total uncertainty with 
respect to a su(2) algebra of each two- level system equals 3/4 — 1/4 = 1/2 and the ratio of 
the uncertainty to the Casimir equals (l/2)/(3/4) = 2/3. Therefore, the strong inequality 
(B7) is not satisfied. More generally, it can be shown using Eq.(B6) that the local algebra 
= ©™ =1 su(d) C su(d n ) gives 

A m i n d 



c n d+V 

therefore the classicality condition (13) does not hold. 
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(B8) 



APPENDIX C: ESTIMATION OF THE NUMBER OF STOCHASTIC REALIZA- 
TIONS, NECESSARY TO CONVERGE THE EXPECTATION VALUES OF THE 
OBSERVABLES IN g TO A PRESCRIBED ABSOLUTE ACCURACY e. 

Given a random variable X with dispersion D x = ^X 2 ^ — ^X^ the number of samplings 
n(e), necessary to estimate the expectation value <^X^ to the absolute accuracy e equals 

n(e) x = ^f. (CI) 

Let us assume that each observable X, e g is measured in an experiment to a prescribed 
accuracy e. The corresponding number of experimental runs is n(e)x 4 - Then 

2 N 



dim{ } V-dim{ } n J2i=l ( X i ) _ ( X i 

2^ n ^ e > x i - — 72 — 



r- e 2 



i=i 

^dim{g} ' 



e 2 



(C2) 



Now consider the computation of expectation values of observables X, e in a state 
evolving according to the Eq. (3), by averaging over stochastic unravelings (26). By Eq.(Cl) 
the number of unravelings necessary to compute the expectation value of Xj to the accuracy 
e is n(e)' x = D' x ./e 2 , where D' x . is the dispersion of the observable in the state p(t). Then 

rdim{g} ' 



E »w> - 1 - 



x 2 - x,- 



v < " e 2 ' " ~ <- 

i=i 



e 2 



(C3) 



where < X >' means statistical average over the unravelings of the quantum expectation 
values obtained in each unraveling (which is the random variable for the purpose of Eq.(Cl)). 
But on the time interval of the simulation (Sec. Ill) 

%)' = (%) , (C4) 



therefore Eqs.(C2),(C3) and (C4) imply 



dim{g} dim{0} 

E = E ( C5 ) 
i=i i=i 
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It follows that 

dim{g} dim{g} 

n{e) st = maxi{n(e)' x .} < n{e)' x . = ^ n(e) Xi < dim{g}max i {n(e) x J 

i=i i=i 

= dim{g}n(e) e:c , (C6) 

where n st (e) is the number of stochastic realizations, necessary to obtain the expectation 
value of each observable Xj G q to an absolute accuracy e, n ex (e) is the number of experi- 
mental runs, necessary to obtain the expectation value of each 5Q to the absolute accuracy 
e. 
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